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The mechanism of lattice diffusion that seems to be best consistent with our col- 
lected experimental knowledge of the metals is that the thermal density fluctuations 
open up voids between the atoms, big enough for one of the adjacent atoms to jump 
into it and that after such a jump the lattice rearranges itself, the net result being 
a change in the mutual arrangement of a few atoms. Such a mechanism is also well 
compatible with a difference in the rate of transport of different atoms in an alloy 
and thus with the Kirkendall effect. In the present paper we shall, however, confine 
ourselves to self-diffusion in pure metals. 

The initiating moment of the process is thus supposed to be the creation of a void 
of the volume w of an atom and the potential energy wu to be stored up in the surround- 
ings of the lattice for this creation is the energy of activation for an elementary process. 
Referred to a mole, the energy of activation is 


E=Nu : (1) 


where JN is the number of atoms of a mole. w is known approximately, and we should 
be able to calculate the energy of activation if we could find the relationship between 
u and w. As the local increase of volume w is supposed to be a result of density 
fluctuations the energy w must also imply the energy of lattice contractions in the 
neighbourhood. 

Now it has been shown in two recent papers (1, 2) that the ratio of the mean in- 
crease of potential energy and the mean increase of volume due to combined density 
fluctuations of both signs, the ratio AU/AV, is a very significant feature of the 
lattice, influencing the change of state on compression, on the thermal expansion of 
solid and liquid and on fusion. The present writer has therefore looked for a relation- 
ship between w/w on one side and AU/AYV on the other. It has not been possible so 
far to derive a relationship mathematically on basis of fluctuation statistics but it 
has been found by trial, that the simple relation 

Emad (2) 
w AV 


leads to very satisfactory agreement with experimental result. A combination of the 
equations (1) and (2) gives 
AU 


= -— Viv, 3 
H=77Nw (3) 
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For w we have to choose the smallest volume that can admit an atom. The opening 
of larger voids will need a higher energy of activation, and as the frequency of a 
certain occurrence of this kind decreases very rapidly with increasing energy of 
activation the lowest value of w will be the decisive one. Thus w is probably the 
volume of the core of the atoms. An apparent value of this volume can be obtained 
from the compressibility formula introduced in an earlier paper (1) 


(Ap +a) (AV +B) =aB (4) 


where Ap and AV are combined variations of pressure and volume, and « and f are 
pressure and volume parameters characteristic of the metal. V and f are suitably 
referred to one mole. With increasing pressure the decrease of volume approaches 
the value -A V =f, and V —f thus represent the volume of the cores of the atoms 
of one mole. We thus have 


Nw=V —f. (5) 
Introduced in equation (3) this gives 
AU 
H=-——(V—-B§). 
ap -B) (6) 


This equation can be used to evaluate the energy of activation H from experimental 
data of various kinds. The formula 


Aes 
AV 2" (" 


deduced in an earlier paper (1) transforms equation (6) into 
3 
ee) (8) 


which enables the calculation of H solely from compressibility data. As such data are 
not yet very accurately known we prefer determining H from equation (6) using 
calorimetric and dilatometric data. 


A suitable measure of the potential energy U is, as shown in an earlier paper (2) 
the total energy minus the kinetic energy expressed by a Debye function, that is 


pif 


3 T 
u-[car—Srrv(z) (9) 


where C' is the molar heat at low pressure, 7’ the absolute temperature, R the gas 


constant, 9 the Debye characteristic temperature determined from specific heat at 
low temperatures, and 


6 
pu 


fi T\3 ( ada 
D(5)=3(5) ( e—1 (10) 
0 
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Fig. 1. Potential energy versus thermal change of volume for tungsten. 


is a Debye function going from 0 to 1 with increasing value of 7'/6. By means of these 
equations and tables for the Debye function, numerical values of U have been ob- 
tained for different temperatures from calorimetric data. 

A plot of U versus the volumes for corresponding temperatures as shown in 
Fig. 1 for tungsten can be used for the determination of AU /AV. As the slope of 
the curve decreases slightly with increasing volume we conclude from equation (6) 
that the energy of activation EH varies slightly with temperature. In order to obtain 
from equation (6) values of # that are comparable with the values obtained from 
diffusion measurements, we therefore have to determine the slope AU/AV from a 
part of the curve corresponding to the range of temperature for these diffusion 
measurements. 

The second factor (V —f) in equation (6) is the volume of the atom cores in a 
mole. It can be supposed to be independent of temperature and thus determinable 
for instance at 0°C. As shown earlier (1) 6 can be replaced by 

eas 


ae wee Vi 1 
p=Vax a Ay” (11) 


where x is the compressibility coefficient. Equation (11) presupposes the use of ra- 
tional units. If, in order to use data from the literature, we take U in calories pro mole, 
V in cm?/mole and x in (kp/cm?)- we have 


(p55 pau 12 
yp 284 Ge (12) 
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Table 1. Data for calculation of the free volume at OC, 


AU / =a Vo x 108 fe 
Vo cm?/mole (kp/em?) + ei 
keal/mole at 0°C 
Li 20.3 13.0 8.8 0.39 
Na 18.0 23.7 15 0.32 
WwW 264 9.6 3.2 0.25 
Cu 60.5 ofp! 0.8 0.19 
Ag 59.6 10.2 1.0 0.16 
Au 73.3 10.2 0.7 0.14 
Al 43 9:9 1.4 0.17 
Fb 38.5 18.3 2.5 0.15 


For a convenient calculation of # from experimental data we finally transform 
equation (6) into the form 


AU AU x 
Seed Borer: 28 a 13 
Bat tec E q ai artad ( ) 
at ¢°C at 0°C 


‘Vos Vo 


which also enables us to refer all changes of volume to the volume V5 at 0°C, as is 
usually done. 

For the evaluation of H by means of equation (13) we need, apart from the com- 
pressibility coefficient x, only diagrams like that of figure 1 for U as a function of 
(V — V,)/ Vo. In the following we shall confine ourselves to eight cubic metals for 
which the activation energies for self-diffusion are determined more or less accurately, 
namely, the body-centered cubic metals Li, Na, and W, and the face-centered Cu, 
Ag, Au, Al, and Pb. For all these metals except Li and W diagrams of U as a function 
of (V — Vo) Vo are given, together with references to the experimental sources, in the 
preceding paper (2). The diagram for W in Fig. 1 is based on tables of specific heat 
at low temperatures by Kelley (3) and at high temperatures by Stull and Sinke (4). 
The thermal expansion data are from Hidnert and Sweeney (5) up to 500°C, and from 
Worthing (6) for higher temperatures. For a diagram for Li we have again used the 
tables of calorimetric data by Kelley (3) and Stull and Sinke (4). The thermal ex- 
pansion data are from the measurements of Pearson (7) and Bernini and Cantoni (8). 

In Table 1 we have collected the data for evaluating the last term in equation (13), 
which according to (12) is 6 / V at 0°C. We can ignore the fact that the compressibility 
coefficients x are mostly determined at room temperature instead of 0°C. The values 
obtained in this way for the free part of the atomic volume, f/V, do not differ too 
much from the values obtained from melting data in the earlier paper (1) for the face- 
centered cubic metals Cu, Ag, Au, Zn, and Pb. For Li and Na, however, the new 
values, 0.39 and 0.32, are much smaller than the values 0.65 and 0.43 obtained from 
melting data. In fact, it was found in the second paper (2) that the melting data give 
too high values for alkali metals, probably because of a structural change in connec- 
tion with fusion. 


Table 2 gives in the first column the mean temperatures for the ranges in which 
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Table 2. Data for calculation of the energy of activation Z. 


e wee [- V V—-B EB B 
Vo ia calc. obs. References 
kcal/mole kcal/mole keal/mole 
Li 100 18.9 0.61 11.5 10 (9) 
Na 50 17.5 0.68 12 10.5 (10) 
WwW 2200 196 0.75 147 142 (11) 
Cu 850 58 0.81 47 48 (12) 
Ag 850 52 0.84 44 46 (12) 
Au 850 61 0.86 52 52 (12) 
Al 500 40.4 0.83 33 33 (13) 
Pb 220 36 0.85 30 28 (14) 


the activation energies are determined. The second column gives the value of the 
first factor of equation (13) at this temperature. The third column gives the second 
factor, and the fourth the calculated value of H. The fifth column gives observed val- 
ues of H from the literature. The values for Na, Cu, Ag, Au, and Pb are from self- 
diffusion measurements, the one for Li from nuclear resonance data, and those for W 
and Al are extrapolated from diffusion measurements on dilute alloys. 

The differences between calculated and observed values for H# are well within the 
limits of the experimental errors. It is hoped that these results will encourage attempts 
to examine whether the gap left in the present paper, a theoretical foundation of 
equation (2), can be filled by a statistical treatment of the fluctuations. 


The author is indebted to Mr. L. E. Larsson for help with numerical work. 


Stockholm, Physics Department, Royal Institute of Technology. 
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